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We analyze a superconductor- ferromagnet (S/F) system with a spiral magnetic structure in the 
ferromagnet F for a weak and strong exchange field. The long-range triplet component (LRTC) 
penetrating into the ferromagnet over a long distance is calculated for both cases. In the dirty limit 
(or weak ferromagnetism) we study the LRTC for conical ferromagnets. Its spatial dependence 
undergoes a qualitative change as a function of the cone angle i?. At small angles ■& the LRTC 
decays in the ferromagnet exponentially in a monotonic way. If the angle l? exceeds a certain value, 
the exponential decay of the LRTC is accompanied by oscillations with a period that depends on 
O ' •&. This oscillatory behaviour leads to a similar dependence of the Josephson critical current in SFS 

»' j , junctions on the thickness of the F layer. In the case of a strong ferromagnet the LRTC decays over 

O ^ the length which is determined by the wave vector of the magnetic spiral and by the exchange field. 
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I. INTRODUCTION 



It is well known that superconductivity and ferromagnetism are antagonistic phenomena (see, for example, the 
reviews QI^E!)- Exchange interaction in ferromagnets results in ordering electron spins in one direction, whereas 
superconducting correlations in conventional superconductors lead to the formation of Cooper pairs with opposite 
spins of electrons. The antagonistic character of ordering in ferromagnets (F) and in superconductors (S) is the 
reason for an essential difference between the proximity effects in S/N and S/F structures (here N denotes a normal 
nonmagnetic metal). In S/N structures the condensate penetrates the normal metal N over a rather long distance 
| which in the dirty limit (tT << 1, r is the elastic scattering time) is equal to 
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6v = y/D/2-KT (l) 

where D — vl/3 is the diffusion coefficient and I — vt is the mean free path. On the other hand, the depth of the 
condensate penetration into a ferromagnet in S/F system is much shorter 



& = y/D/h (2) 

if the exchange energy h is larger than the temperature T (usually h >> T). The formula 10) is valid for a short mean 
free path ( hr << 1). If the exchange field h is strong enough (hr >> 1), the condensate penetrates the ferromagnet 
over a distance of the order of I (if tT < 1) and oscillates with the period ~ v/h 0,O|- 

Note that a short length of the condensate penetration is related to the fact that Cooper pairs in a conventional 
superconductor are formed by two electrons with opposite spins. In the case of a homogeneous magnetization the 
wave function f(t — t') of Cooper pairs penetrating into the ferromagnet consists of two parts 



P., h(t- 1') ~ (Vt(W) - Mt')^(t)) (3) 

fo(t - ~ WtWiV) + Mt')^(t)) (4) 

The first function describes the singlet component. It differs from zero both in the whole superconducting region 
and in the ferromagnet over the length £p. The second function fo(t — t') describes the triplet component with zero 
projection of the magnetic moment of a Cooper pair on the z-axis (the magnetization vector M is oriented along the 
z-axis). It is not zero only in the vicinity of the S/F interface, over a distance ~ £s = yj D /2nT c in the superconductor 
and over the distance £p in the ferromagnet. This function is an odd function of the difference (t — t') and therefore 
is equal to zero at t = t'. This means that in the Matsubara representation fo(u>) is an odd function of u>, whereas 
the function /^(w) is an even function of lo. 
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All the statements above concern only the case of a homogeneous magnetization in the F region. The situation 
changes qualitatively if the magnetization is not homogeneous, for example, if one has a spiral magnetic structure in 
the ferromagnet. In this case not only the singlet and triplet component with zero projection of magnetic moment, 
but also a triplet component with a nonzero projection of magnetic moment arises in the system. This type of the 
triplet component means that Cooper pairs appear in the system that are described by the condensate function 
f tr (t - t') ~ (^tW-0t(*')) or ftr(t - t') ~ (tp L (t)?p L (t')). The triplet pairing is well known in superfluid He 3 d Q 
or in Sr2RuO^ |10| and is believed to be realized in compounds with heavy fermions The condensate function 
ftr(w) in these materials is an odd function of momentum p and an even function of the Matsubara frequency u>. The 
order parameter A is a sum of ftr{ui) over positive and negative u> and corresponds to a triplet pairing. This triplet 
component is suppressed by impurity scattering [T^ . 

In S/F structures the impurity scattering is rather strong as ferromagnetic films used in these structures are thin 
(typically the thickness of the F films d is about 20-100 A) and the elastic scattering at least at the F surface is 
strong. Therefore the conventional triplet component would be strongly suppressed. However there is a special type 
of the triplet component which can survive a strong impurity scattering. This triplet component is described by a 
condensate wave function that is even in momentum p and odd in frequency lu. This type of the condensate 

was first suggested by Berezinskii in 1975 in attempt to describe the pairing mechanism in superfluid He 3 [13J. It 
turned out however that in reality the condensate function in He 3 is an even function of u> and an odd function of 
p. Later a possibility to realize the odd (in frequency) triplet superconductivity in solids was discussed for various 
models in Refs. jl4L fl5l Il6| ■ 

The odd in frequency u> and even in momentum p triplet component in S/F structures differs from that discussed 
in the preceding paragraph. It coexists with the singlet component and the order parameter A is determined only 
by the ordinary (BCS) singlet component / S n ff z(w) even in uj. The triplet component with nonzero projection of the 
magnetic moment of Cooper pairs arises as the result of action of a rotating exchange field on electron spins. This 
type of the condensate penetrates the ferromagnet over a long distance of the order of £jv (see Eq.Q) provided the 
period of the magnetization rotation exceeds £n- 

A S/F structure with an inhomogeneous magnetization has been studied for the first time in Rcf. 6J. The authors 
considered a S/F structure with a Bloch-type domain wall at the S/F interface in the limit of a short mean free 
path (Jit < 1). In the domain wall of the thickness w the magnetization vector was supposed to have the form 
M =Mo{0, sinQx, cos Qx}, where the x-axis is normal to the S/F interface. Outside the domain wall the magnetization 
was constant: M = Mo{0,sin(5w,cosQw}. The condensate function f tr (ui,x) ~ {ipa{t)' l Pa{t , )}u> was found from the 
linearized Usadel equation. It was established that this triplet component odd in frequency and even in momentum 
penetrates the magnetic domain wall over the length 



t Q = {Q 2 + 2n\u\/D}- 1 / 2 (5) 

It spreads outside the domain wall over distances of the order £jv. At x > w the vector of the magnetization is fixed 
so that the first term in Eq.@J Q 2 should be dropped. This triplet component may be called the long-range triplet 
component (LRTC). The LRTC may cause a significant change in the conductance of the Andreev interferometer 
consisting of ferromagnetic wires and a superconducting loop. As was shown in Ref.[6(, the conductance variation 
decreases with increasing temperature in a monotonic way. 

Somewhat later the same problem was considered in the paper |17| . In that publication a more complicated situation 
was discussed, namely, the case when the width of the domain wall is short compared to the mean free path. In order 
to find the condensate function (quasiclassical Green's function) f tr (oj,x) in this case, one needs to solve a more 
general Eilenberger equation taking into account a non-homogeneous magnetization. Unfortunately, the authors of 
Ref.^3 did not manage to solve the Eilenberger equation and therefore restricted themselves with a rough estimation 
of the amplitude of ft r (ui,x). So, the problem of calculation of the odd triplet condensate function in the ballistic 
regime remained unsolved. 

In the present paper we continue studying behavior of the triplet condensate considering new situations. To be 
specific, we study the LRTC in S/F junctions with a spiral magnetic structure in different limits including the quasi- 
ballistic one (i.e. hr > 1). This spiral structure may be both an intrinsic property of a ferromagnet (for example, a 
helicoidal structure; see [l8j) or may just serve as a rough model for magnetic domains. 

We consider a system with a conical ferromagnet, that is, we assume that the magnetization in F rotates in {y, z}- 
pla ne and has a constant component along the x-axis. This type of spin structures is realized, for instance, in Ho 
|19| and what we discuss now is a generalization of the problem considered in 0, where the ^-component of the 
magnetization in F was assumed to be zero. Our study of the proximity effect in such a S/F system with a spiral 
magnetic structure is motivated also by the recent experiment performed on a Al/Ho structure p0| . 

First, solving the Usadel equation, we find the condensate function in the dirty limit. It turns out that the spatial 
variation of the LRTC has a nontrivial dependence on the cone angle ■& (see Fig.l). If ■& is small, the LRTC decays 
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with x exponentially over a distance of the order of £q, but at angles > sin - (1/3) « 19°, the LRTC decays in a 
non-monotonic way. It oscillates with a period depending on i?. These oscillations lead to oscillations of the critical 
Josephson current in SFS junction as a function of the thickness of the F film 2L if the thickness 2L is essentially 
greater than £p. In this case the Josephson coupling is only due to the LRTC. 

We continue the investigation studying in the third section the LRTC in a S/F structure with a spiral structure in 
the limit when the condition hr >> 1 is valid and the Eilcnberger equation must be solved. We analyze peculiarities 
of the LRTC in this case. Surprisingly, this case has not been investigated previously, although it may correspond to 
a real situation. In many ferromagnets the exchange energy h is very large so that the product hr can be arbitrary. In 
both sections we make an assumption that the proximity effect is weak, i.e. the amplitude of the condensate function 
in the ferromagnet / is small. This assumption is presumably valid in most cases because there is a strong reflection 
of electrons at the S/F interface due to a considerable mismatch in electronic parameters between the ferromagnet 
and superconductor |29f . In conclusion we summarize the obtained results and discuss a possibility of experimental 
observations of the predicted behavior of the LRTC. 

Note that in Refs. 0,0] and in the present paper a Bloch-type spiral structure is analyzed. The rotation axis was 
assumed to be perpendicular to the S/F interface so that the condensate decays in the ferromagnet in the direction 
parallel to this axis (see Fig.l). In this case the LRTC arises in the system. 

At the same time, one can imagine another, Neel-like type of a spiral structure with the magnetization vector M 
that rotates in the plane of the ferromagnetic film and does not vary in the perpendicular direction. A solution for 
this type of the spiral structure has been found in Ref. [23 • The authors considered the case of a thin F film and did 
not study the decay length in the direction normal to the film. 

The same problem for a thick (infinite) ferromagnetic film was analyzed by Champel and Eschrig in a recent paper 
[23| . These authors assumed that the magnetization vector M lies in the plane of the ferromagnetic film and rotates 
in this plane being constant in the perpendicular direction. Surprisingly, they found that there was no LRTC in this 
case. 

As was shown later [24|], the absence of LRTC is specific for this homogeneous spiral magnetic structure. In contrast, 
the long-range triplet component arises if the Neel-type spiral structure is non-homogeneous; for example, if the F film 
consists of magnetic domains separated by the Neel walls. In this case the LRTC appears in domain walls and decays 
in domains over a long length of the order £at. Another case where the LRTC arises is a spin-active S/F interface. 
Such type of the interface leads to mixing singlet and triplet component and the triplet component may penetrate 
even into a half-metal ferromagnet when the conduction band consists of electrons with only one direction of spins 

II. S/F STRUCTURE WITH A ROTATING MAGNETIZATION. DIRTY LIMIT 

We consider a system shown schematically in Fig. 1. A ferromagnetic film F is attached to a superconductor S. 
The magnetization M in the ferromagnet is assumed to rotate in space and the vector M has the form 

M =Af {sin cost? sin Qa;, cost? cos Qx} (6) 

This form of the magnetization vector implies that the projection of the M vector onto the x-axis Mo sin $ is constant 
and the projection onto the {y, z}-plane Mo cos?? rotates in space with the wave vector Q. We also assume that the 
condition 

hr « 1 (7) 

is fulfilled, that is, cither the exchange energy h is not large or the collision frequency t _1 is high enough. In addition, 
we assume that the proximity effect is weak, i.e., the condensate function / is small. The smallness of the / function 
means the presence of a barrier at the S/F interface or a big mismatch in electronic parameters of the superconductor 
and ferromagnet (the Fermi momenta in F and S differ greatly). In this case one can linearize the Usadel equation 
and represent it in the form 0, Q 

d 2 f/x 2 ~2klf + ik 2 h sgnu;{cos$([a 3 j] + cosa{x)+T 3 [a 2 , f] + sma(x)) + sini9f 3 [a u f]} = 0, (8) 

where u> = irT(2n + 1), k 2 , = \uj\/D, k 2 = h/D, the brackets [173,/]+ and [03,/] denote anticommutator and 
commutator. The quasiclassical condensate Green's function / is a 4 x 4 matrix in the particle-hole and spin space. 
The Pauli matrices fi and &i operate in the particle-hole and spin space, respectively. 
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FIG. 1: Schematic view of the system under consideration, a) S/F structure with a conical ferromagnet. Magnetization in the 
ferromegnet F rotates in the {y,z} plane and has a constant projection on the x-axis. b) Josephson junction with a conical 
ferromagnet. 



The Green function in the bulk of the superconductor is 



fs = 7~2 &3 fs 7 



(9) 



with fs — A/vw 



A 2 . This function describes an ordinary singlet condensate in conventional superconductors. 
Assumingthis solution to be valid up to the interface, we use the "rigid" boundary condition at the SF interface 

( X = o) mm 



df/dx = -fs/lB, 



(10) 



where jb — Rb&, & is the conductivity of the ferromagnet and Rb is the interface resistance per unit area. This 
boundary condition is valid provided the ratio £h/jB — p£,h/RB is large, that is, the interface resistance should be 
larger than the resistance of the F wire of length 

What we have to do is to solve Eq.© in the ferromagnetic region (x > 0) with the boundary conditions Q OJl at 
x = 0. Eq.(JSJl is a linear differential equation with coefficients depending on the coordinate x. This dependence can 
be excluded if we make an unitary transformation determined by a matrix U and introduce a new matrix /„ 



/ = Uf n U+, 



(11) 



where the matrix U = exp(iT3<7iQx/2) describes a rotation in the spin and particle-hole space. Substituting the 
expression (|ll|l into Eq.(|SJ), we get the equation for the new matrix /„ 



~d 2 f/dx 2 + (Q 2 /2 + 2kl)f + (Q 2 /2)(a 1 fa 1 )^tQf 3 [a 1 ,df/dx} + -ik 2 h sgnu J {cosd [a 3 , f] + + sin tff 3 [a 1: f]} = 0, (12) 
For brevity we dropped the subindex n. The boundary condition (|10|) acquires the form 



df/dx + t(Q/2)fa [a 1 ,f] + = -fa h B > (13) 

As concerns the matrix fs, it does not change (it is invariant with respect to the transformation l|ll[l 'l. In order to 
solve Eq. i|12|) with the boundary condition (|13|) . we represent the matrix f(x) as an expansion in Pauli matrices 



f{x) =^2&iFi{x) 

i 



(14) 
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where i = 0,1,2,3 and <7o is the unit matrix. The functions Fi(x) are matrices in the particle-hole space. The 
spatial dependence of the functions Fi(x) is determined by the exponential functions Fi(x) ~ Aiexp(— kx), where 
the eigenvalues k are determined from the determinant of Ea. (|12|l . Putting the determinant to zero, we obtain the 
equation for eigenvalues k (see Appendix A) 

(k 2 - Q 2 - 2k 2 ) 2 [(k 2 - 2k 2 ) 2 + 44 sin 2 1?] + 4(Qk) 2 [(k 2 - 2k 2 ) 2 + 4«# sin 2 1?] + 4«# cos 2 i?(k 2 - Q 2 - 2k 2 ) (n 2 - 2k 2 ) = 

This equation has 4 pairs of roots. However, four of them correspond to solutions growing in the ferromagnet. As 
we are interested only in decaying solutions (Rck > 0), we keep four proper roots. 

In order to simplify the calculations, we consider a limiting case of a large exchange field assuming that 

k 2 »{Q 2 ,k 2 J (16) 
In this limit two eigenvalues of the wave vectors k are large and equal to 

k± = (1 ± isgnu)Kh (17) 

These roots describe a rapid decay of the condensate in the ferromagnet over an "exchange length" k^ 1 of the order £p 
(see Eq.<0). Only these eigenvalues appear in the case of a homogeneous magnetization. These large eigenvalues have 
both real and imaginary parts and they are equal to each other. This means that corresponding eigenfunctions decay 
and oscillate in space on the same scale. The oscillation of the function F3 (singlet component) lead to oscillations of 
the critical temperature T c and Josephson current I c in SFS structures with varying the F film thickness d 0, S Q • 
In addition to these large values K±, there are two other solutions K a ^ that describe a long-range penetration of the 
triplet component. The singlet component also contains a part decaying slowly but, as we will see, its amplitude is 
small in comparison with the amplitude of the triplet component. These eigenvalues are small: n 2 a b sa max{Q 2 , k 2 }. 
One can find exact expressions for K a ,b, but in order to make results more transparent, we represent solutions in the 
limit 

Q 2 » *l (18) 
In this case Ea. 1)15(1 is reduced to the following quadratic equation for the eigenvalues ( a ,b = k 2 b /Q 2 

C 2 - C(l - 3sin 2 tf) +sin 2 tf = (19) 

with the roots 



Ca.b = 3/2 • [1/3 - sin 2 ± cos 1?^/ 1/9 -sin 2 ??] (20) 

Now, let us consider these eigenvalues as a function of the angle 1?. If the magnetization in the ferromagnet lies in 
{y, z} plane (i.e., z? = 0), we obtain only one root: £ a = 1 (the other root should be dropped as it corresponds to a 
solution with the zero amplitude). This means that the triplet component penetrates the ferromagnet over the long- 
range distance (^lr ~ Q^ 1 . If the value of Q 2 is comparable with or less than k 2 , then the characteristic penetration 
length of the triplet component is £la = 1/ \JQ 2 + (compare with Ref. @). As follows from Ea. f2T)jl . at small 1? 
the first eigenvalue equals ( a ~ 1 — 4 sin 2 •& ps 1 — 4?9 2 , whereas the second root is small and equal to ~ sin 2 •& m ?? 2 . 
Therefore in this case the exponential decay of the triplet component is slow: ft r ~ exp{—Qx\'&\). 

Both the eigenfunctions decays exponentially in a monotonic way over a length much longer than £p . The situation 
changes if sin $ exceeds the value 1/3 (8 > sin -1 (1/3)). In this case the LRTC oscillates and decays in a non- monotonic 
way. The period of the oscillations ~ 1/Q is much longer than £p. If the magnetization vector is oriented almost 
along the cc-axis (sintf — > 1, cos?? — > 0), the period of oscillations ~ Q^ 1 is much shorter than the decay length of the 
LRTC which is equal to V2/(Qcostf) » 1/Q (see Fig.2). 

The amplitudes of the components F(x)i are found from the boundary conditions l|13fl . One can show that the 
matrices F(x)i can be represented in the form (see Appendix) 
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F(x)o,3 = t 2 2J A a,3k exp(-K fc x) (21) 



T3F(x)x t 2 = f 2 ^ Al - 2k exp(-K k x) (22) 

where the summation is carried out over all eigenvalues: k = ±; a, b. In the approximation, Eq. I|16|). the amplitudes 
are 

T A Q± / cosd = ±iA 2 ±/ sintf = A 3 ± = (/ s /2 7flK± ) (23) 



-itwdA 2a ^ b = A Qafi = (C a , 6 - l)/(2iVW0^ia,& (24) 
The amplitude of the long-range triplet component is equal to 

, fs CaCOSZ? 

A la = sgnw =- (25) 

7B>«h (Ca + 1)(VU-V5) 

The expression for An, is obtained from Eq. l|25l) by permutation a <=^ b. These formulas are valid if both sin$ and 
cost? are not too small: {cos^, sind} » Q/kh. 

As we mentioned the amplitude A3 corresponds to the singlet component and the amplitude Aq describes the 
triplet component with zero projection of the magnetic moment on the z— axis. In a homogenous case (M is constant 
in the ferromagnet) this component penetrates the ferromagnet over a short length of the order £p. In the case of 
rotating magnetization it penetrates the ferromagnetic region over a long distance of the order of min{l/Q, 1/ac w }. 
The amplitude A 2 arises only in the case of a tilted magnetization (d 7^ 0). One can see from Eas. (|23l25[) that 
in the considered limits the amplitudes A$±, Aq±, A 2 ±, -Aoa,&, A 2a ,b and A\ a ^ are comparable with each other (other 
amplitudes are small). Therefore the LRTC with nonzero projection on the z— axis is comparable with the magnitude 
of the singlet component at the S/F interface (^3+ + A3-) ~ fs/l K h- The spatial dependence of the LRTC amplitude 
is given by the expression 



where /lr{ x ) is defined in this way: Fi(x) — ififLn(x), x » £p. 

As it should be, this function is an odd function of ui. In Fig. 2 we plot the spatial dependence of Re(/ifl(a;) for 
some values of sin$ and compare it with the spatial variation of the singlet component. 



/si 1 

h( x ) = o — [ — exp(-K+x) H exp(-K-x)} (27) 

27s «+ «- 

It is seen that the LRTC decays over distances longer than the singlet component. At small sin$ the LRTC decays 
monotonously, but with increasing oscillations of the LRTC arise. The characteristic scale of the LRTC decay in 
the considered case (the condition (|18|l is fulfilled) is Q" 1 . In the opposite limit of a spiral with a small Q (Q << nj) 
the roots K a ,b changes: n a ,b ~ V / 2'*w(l ± iQ sin '&/\2k U i)- This means that the LRTC decreases exponentially in the 
ferromagnet on the length of the order £jv and oscillates with the period (Q sin-!?) -1 . In this case the decay length is 
shorter than the oscillation period and therefore there should be no oscillations in observable quantities. In a general 
case the characteristic length l tr of the LRTC decay is k r ~ min{l/Q, £jv}. 

In the next Section we discuss a possibility to observe the unusual behavior of the LRTC and demonstrate that 
such an observation is realistic. 
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III. JOSEPHSON EFFECT IN SFS JUNCTION 



In this Section we consider a Josephson SFS junction with the same spiral magnetic structure as in the preceding 
Section. This junction is shown schematically in Fig. lb: a ferromagnetic layer with the spiral structure connects two 
superconductors where the phases of the order parameter are equal to ±<p/2. We assume again a weak proximity 
effect and consider an important case where the Josephson coupling between the superconductors is only due the 
LRTC only. This means that the condition 



n h L » 1 (28) 

should be fulfilled. In this case the singlet component decays fast near the S/F interfaces and the Josephson coupling 
is provided by an overlap of the LRTC. 

In order to find the Josephson current Ij, we need to solve Eq. l|12|l with boundary conditions at x = ±L. The 
boundary conditions coincide with Eq. \i'A\ . where in this case the matrix fs at x = ±L is replaced by 



fs = /s(r 2 cos(^/2) ± h sin^/2)) • a 3 (29) 

The matrix f(x) is again represented in the form of an expansion in the Pauli matrices (see Ea. l|14|) '). but the 
"coefficients" F(x)i in this expansion acquire a somewhat more complicated form. These matrices (in the particle- 
hole space) have the structure (see Appendix B) 



F(x)o,3 = h ^ A ^,3k cosh(K k x) + ti ^ A ,3k sinh(K fe a;) (30) 
k k 

The matrix fsF(x)2 has the same structure. As concerns the matrix T3-F(x)i, it has a similar structure but with 
another spatial dependence 



hF{x)i — f 2 ^2 A lk sinh(K fc a;) + n ^ Aim cosh(K k x) (31) 



k 



The coefficients Ai k ,Aik can be found from a solution for Ea. l|12|l in a way similar to that in the preceding Section. 

The Josephson current (per unit square) is calculated using a general formula for the condensate current in the 
dirty case (see, for example, [71127)1 
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I j = i(nT/Ap) Tr{r 3 ■ a fdf/x} (32) 

OJ 

where p is the resistivity of the F metal. Making the transformation , one can rewrite this formula in terms of 
the new function /„ 

I j = (inT/Ap) Y, Tr{(t 3 -a )Kdf/x + iT 3 (Q/2) [/>!]+)} (33) 
We dropped again the index "n". After simple but somewhat cumbersome calculations we obtain for Ij 

Ij = I c am<p, I c = nT(a/2)\Q\ Y^{A la A la {\ - ^7^ )/^} (34) 

' 4sm v 

Lj.a 

where the summation over a means that a = a and a = b. 

The amplitudes Ai a ,Ai a can be found as before. We find (see Appendix B) 

, Ka i{A 3 - -A 3+ )cosd 

Ala - (C„ + l) S inh0 a M [6b) 

. _ 2Ca i{A 3 --A 3+ ) COST) 

• Ala - (Ca + lJcOSh^ M [ ' 

where M — VCa/ tanh# a — y^/tanh^b, Ai = \/~C2ta, i ahO a — \fQ, tanh Of,, 6 a ^ = K a ,bL. The coefficients Aib,Aib are 
given by the same formula with replacement a — > b. Under the condition l|28[) the coefficients for the singlet component 
A 3 ± arc described by formulas similar to Eq. (|23|) 

A 3± = -^— cos^/2); A 3± = — ^— sin(<^/2) (37) 
2-f B K± 2j b k± 

Eqs. H34I37J) determine the Josephson critical for the SFS junction under consideration. In the approximation l|18|l 
we can perform the summation over u) and obtain for the critical current 



2 AQ\ A W A 

m Atanh( — , 



Ic= „,.,_„,\ 2 Atanh( — )J C (38) 



Here the coefficient ^b^k is assumed to be large. Only in this case the Usadel equation may be linearized. However 
the obtained results are valid qualitatively in the case when this factor is of the order 1. The quantity J c depends 
only on the angle and the product QL. It is equal to 

_ cos 2 ^ ^ C, a 2 {t a ~l)\ 1 , , , 

Here the roots £ a .b are given by Ea. (|20|l . Since at some angles the roots have an imaginary part, one can expect 
that the normalized critical current changes its dependence on L with varying angle z9. First we demonstrate this 
analytically considering a limiting case. We assume that the overlap of the condensate induced by each superconductor 
is weak. This means that \6 a \ >> 1 and therefore sinh26 a j, w (1/2) exp(20Q,,b), tanh# Q sw 1. If in addition the angle 
$ is close to 7r/2,i.e. cosi? << 1. In this limit we obtain 

cos $ exp (-V^Q£cosi?)sin(2QZ) (40) 
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FIG. 3: The dependence of the normalized critical Josephson current J c on QL for various S = sini}. 

Therefore the normalized critical current J C (L) as a function of L undergoes many oscillations with the period (tt/Q) 
on the long decay length (y^Qcosi?) -1 . It turns to zero at cosz? — > 0, but the decay length becomes infinite. We 
remind that there is a lower limit on cos$. It was assumed that cosd is larger than the small ratio Q/nh- The 
maximum of J c is achieved at cos <d w 1/V2QL if QL » 1. 

In Fig. 3 we plot the dependence of the critical current I c on the length L for different projection of the magnetization 
on the x-axis. In accordance with the analysis above, it is seen that at the angle determined by sin?? = 1/3 the decay 
of I c with increasing L is accompanied by oscillations. The period of these oscillations is of the order Q _1 and depends 
on the angle d. These oscillations are caused by oscillations of the LRTC. 

IV. STRONG FERROMAGNET (QUASI-BALLISTIC CASE) 



In this Section we consider the opposite limit, i.e. the condition 

hr » 1, (41) 

is assumed to be fulfilled. This case may be realized either in a weak ferromagnet with a large mean free path or in 
a strong ferromagnet with a large exchange energy h. For example, this condition is fulfilled for h « leV if the mean 
free path is longer than ~ 50 A (we take the Fermi velocity v = vft, ~ vfi ~ 10 8 cm/s). For simplicity we assume that 
the M vector lies in the {y, z}-plane (that is, # = 0). In order to find the condensate function in the ferromagnet, we 
have to use the more general Eilenberger equation. We assume again that the proximity effect is weak. The linearized 
Eilenberger equation for the condensate matrix function / reads [2£| (see Appendix A) 

Hl{f 3 df/dx + i(Q/2) fa,/] + } + 2Mt/ - ^rsgnu, [a 3 ,f] + = </) - /, (42) 

where /i = p x /p, Px is the projection of the momentum vector p on the cc-axis, I = tv is the mean free path, the angle 
brackets means the angle averaging. The boundary condition is |2l| 

(sgnw)f 3 a = (sgn^)i M / s J a;= o (43) 

where a(/i) = [f(fi) — /(— /x)]/2 is the antisymmetric part of the condensate function and = T([i)/4, T(p) is the 
transmission coefficient of the S/F interface. In order to find the solutions of Ea. (|42|) . we represent the matrix / as a 
sum of antisymmetric and symmetric functions 



f(/j.) = a(n) + s(fi) 



(44) 
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where s(ji) = [/(/x) + f(—fx)]/2. The matrices a{fi) and s(/x) are represented again as a series in the spin matrices &i 



E- 



(45) 



where i = 0,1,3. The coefhcients 

Ct^i <§2 3X6) elS before, matrices in the particle-hole space. We substitute these 
expansions into Ea. H42|) and single out the symmetric and antisymmetric parts. After some algebra we obtain for the 
diagonal matrix elements in the spin space (s)n(22) — sq i §3 = s± (see Appendix C) 



s±(x) = ±T2t fl f s exp(-K±x/ln) (46) 

where K± = =F icx-h, a u = 1 + 2\lj\t, a/i = 2/irsgnoj. The singlet (S3) component and the triplet (sq) component 
with zero projection of the magnetic moment of Cooper pairs are related to s±(x) : 83,0(2;) = (511(2;) =F s 22(2:))/ 2- As 
is seen from Eq.l|46|l. these components oscillate in space with a short period of the order of v/h and decay over the 
mean free path / (if tT << 1). When deriving the expression for ,§3,0(2;), we assumed that the condition 



Ql/\a h \ = Qv/h « 1 (47) 

is satisfied, but the relation between the period Q^ 1 of the spiral and the mean free path may be arbitrary. 

Let us turn to the more interesting LRTC Si. An equation that describes the LRTC differs considerably from the 
one for the matrices s±(x). Characteristic wave vectors are much smaller in this case. For the Fourier transform 
Si(fc) = / dxsi (x) exp(ifcir) this equation has the form 



Wl + {Q 2 (a u /a h f + /c 2 )(^) 2 ]si(fc) + {l^) 2 (Qk)T 3 [so + i(a u /a h )s 3 ] = <*u(si) (48) 

The matrices .§0,3(2;) are given in the first approximation by Ea. (|46l) . Ea. H48|) can be solved in a general case, but 
we are interested in the behavior of the LRTC at distances much longer than the mean free path I. This means that 
one has to find the matrix si(fc) for small k: k << I. We find (see Appendix C). 



^) = ~^^\)fsin (49) 



where (fJ- 2 t^) = d/ifi 2 t^ and Kq = 2\u>\/D + (Q/cth) 2 - Performing the inverse Fourier transformation, we find the 
spatial dependence of the LRTC 



si(x) = -3 — — (i^t^) fs exp(-xK Q )ifi 7 (50) 

OthJ<Q 

We took into account that a 2 ~ 1. 

This is the main result of this Section. The spatial dependence of the singlet component £3(21) = (s+ — s_)/2 
is given by Eq. (|46f) Comparing Eas. (|50(l and l|46|l . we see that the amplitude of the LRTC is comparable with the 
amplitude of the singlet component at the S/F interface. Indeed, if \/2itT / D « Q/cth, then Kq ~ Q/cth and the 
coefficient (Q/a^i^Q) in E!a. (|5U|I is of the order of 1. 

In Fig. 4 we plot the spatial dependence of the LRTC and the singlet component. One can see that the singlet 
component oscillates fast with the period of the order v/h = l/hr (we take the magnitude of hr equal to 5) and 
decays over the mean free path I. The LRTC decays smoothly over a length ~ (Q//ir) _1 (we assumed that (Q/hr) 2 > 
2itT/D). 



V. CONCLUSIONS 



We considered the odd triplet component of the superconducting condensate in S /F systems with a spiral magnetic 
structure. The axis of the spiral is assumed to be perpendicular to the S/F interface (the Bloch-like spiral structure). 
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Si S3 




a) b) 



FIG. 4: Spatial dependence of the real part of the LRTC(a) and singlet (b) component in quasiballistic case. The singlet 
component shown for fi = 1 does not depend on Q. Both component are normalized to the quantity (/x 3 t M )/s and {t^fs ), 
respectively. The parameter hv/Q is choosen equal to 5. 



We analyzed both dirty (hr < 1) and clean (hr > 1) limits. These limits correspond in practice to the cases of weak 
and strong ferromagnets. 

In the diffusive limit we studied the case of a conical ferromagnet when the magnetization vector M has the 
constant projection M siwd on the x-axis and rotates around this axis. The condensate amplitude in the ferromagnet 
was assumed to be small compared to its amplitude in the superconductor (a weak proximity effect). In addition, we 
assumed that the exchange energy h was larger than such energy scales as T, DQ 2 (dirty limit), Qv (clean limit). In 
this case the penetration length of the singlet component is much less than that of the LRTC. The singlet component 
penetrates the ferromagnet F over a distance of the order \J D/h in the dirty limit and over the mean free path I in 
the nearly clean limit (if tT << 1). In the ballistic case (tT >> 1) the singlet component decays over distance of 
order v/T. In the nearly clean or ballistic case the singlet component oscillates fast with the period v/Q. The LRTC 
decays over a length of the order min{Q _1 , £at}. 

In conical ferromagnets the LRTC has an interesting nontrivial dependence on the cone angle At small $ the 
LRTC decays exponentially in a monotonic way over the length 1/Q (if Q >> C^ 1 ), but at > sin _1 (l/3) « 19° 
the exponential decay of the LRTC is accompanied by oscillations. The period of these oscillations depends on $ so 
that at i? — ► 7t/2 the period of oscillations is much smaller than the decay length. The amplitude of the LRTC is 
comparable with the amplitude of the singlet component at the S/F interface. The latter amplitude is determined by 
the S/F interface transmittance and decreases with increasing h. 

In the dirty limit we calculated also the critical Josephson current I c for a S/F/S junction with a conical ferromagnet 
F. It was assumed that the thickness of the F layer 2L is much larger than £p. Therefore the Josephson coupling 
is only due to an overlap of the LRTCs whereas the overlap of the singlet components induced by superconductors 
is negligible. The dependence of I c on the angle $ is determined by the LRTC: at small $ the critical current I c 
decreases with L monotonously, but with increasing $ the decay of the function I C {L) is accompanied by oscillations. 
Therefore measurements of the Josephson critical current in SFS junctions with a conical ferromagnet may provide 
useful information about the LRTC. 

Note that the triplet component may also exist in magnetic superconductors with a spiral magnetic structure (see 
Ref.Q and refer ences therein). However in magnetic superconductors the triplet component coexists with the singlet 
one and, contrary to our case, can not be separated from the singlet superconductivity. 

We also studied the LRTC in the limit hr > 1 for -d = 0. In this case the singlet component decays over a length 
of the order of the mean free path I (if tT << 1). Its amplitude at the S/F interface is determined by the S/F 
interface transmittance and does not depend on h. The LRTC penetrates the ferromagnet over a length of the order 
of (Q/hr)^ 1 (if Q > (}it)/£n). The decay length of the LRTC is longer than the decay length of the singlet component 
I provided the condition h > Qv is valid. 

Note that we neglected the spin-orbit interaction. The latter restricts the penetration length of the LRTC by the 
value of the order of D/8t s . . where r s o . is the spin-orbit relaxation time 

We would like to thank SFB 491 for a financial support. 
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VI. APPENDIX A. BASIC EQUATIONS; DIRTY CASE 



The Eilenberger equation in a stationary case for the system under consideration has the form 



dg 



dx 



1 



(j,v- ie[f 3 a ,g]- A<T 3 ,g - ih{cosi)([T 3 a3,g] cosa(x) + [a 2 ,g] siaa(x)) + sintf [<ti, g])} + — [(g), g] = 0, (51 



2r 



where fx = p x /p, v is the Fermi velocity, g is a 4 x 4 matrix of the retarded (advanced) Green's functions. The order 
parameter A is not zero only in the superconductor and the exchange field h differs from zero only in the ferromagnet. 
If the phase in the superconductor is chosen to be zero, then A = if^A. The angle brackets mean the angle averaging. 
In the Matsubara representation the energy e is replaced by w: e — > iui. If the condition hr << 1 is fulfilled, the 
antisymmetric in momentum p part of g(p) is small and can be expressed through the symmetric part of g. For the 
symmetric part g one can obtain the Usadel equation (see, for example, 0) that in the Matsubara representation 
reads 

Dd(gdg)/dx 2 + uo [f 3 <T ,<?] _ Aer 3 ,g + i/i{cos$([T 3 <7 3 , g] cosa(x) + [62,5] sina(x)) + sin$ [cti,*?])} = 0, (52) 

where D — vl/3 is the diffusion constant. 

In the case of a weak proximity effect one can linearize Eas. <|51|) and (|52|l . For example, in order to obtain the 
linearized Usadel equation for a small condensate function f(x) in the ferromagnet, we represent g(x) in the form 

g(x) = sgnwf 3 o- + f[x), (53) 

where the first term is the matrix quasiclassical Green's function of a normal metal. Linearizing Eq. l|51|) with respect 
to f(x), we come to Eq.@. 

In order to find solutions for the matrix f(x), we substitute the representation of this matrix f(x) in the form 
of (|14fl and 121I22|I into Eq. (|12|) . As a result we obtain on the left-hand side of this equation a sum of four terms 
proportional to the matrices <7;. Coefficients at each matrix &i are matrices in the particle-hole space. The sum of 
these four terms equals zero. Therefore we obtain four equations for these coefficients at &i, where i = 0, 1, 2, 3: 

(-k 2 + Q 2 + 2kI)F {k) + 2iQ«(faA(«)) - 2icostiK 2 h F 3 {K) = 0, (54) 
2iQkF {k) + {-k 2 + Q 2 + 2kI){t 3 F 1 (k)) = 0, (55) 
(- K 2 + 2kI)(t 3 F 2 {k)) -2sm#K 2 h F 3 (K) = 0, (56) 

- 2icosi9/c 2 [ F (/t) + 2sin^(f 3 F 2 (K)) + (-k 2 + 2k 2 )F 3 (k) = 0, (57) 

This system of equations has a nonzero solution if the determinant of the system is zero. Thus we come to Eci. <|15[) 
for four eigenvalues Kk- In order to determine the matrices Fi, we have to use the boundary conditions (|13(l . Since 
the matrix fs in the particle-hole space contains only the matrix T2, the matrices Fq_ 3 (k) and t 3 Fi^(k) also are 
proportional to the matrix T2. Therefore one can write the expansion l|21l) - (|22l) . The amplitude of each mode are 
Ai y k, where the first index i is related to the spin space and the second k mean the eigenvalues of the wave vectors Kk 
(k = ±,a, b). From Eas. (|54l57|l one can determine relations between amplitudes Ai : k(nk) for each mode. As follows 
from Eqs. I|54I57|I . the coefficients Ai(nk) are connected with each other via Eas. l|23|l - H25|) provided the conditions illtj|) 
and (fTHfl are fulfilled. The matrices (t 3 Fi(k±)) w (2iQ/ k±)F (k±) and F 3 (K a ,b) ~ —{Q/Kh) 2 {z a ,b/ sint?) 2 F (k±) are 
small compared to other matrices. Using these relations and substituting the representation (|14|l into the boundary 
condition (|13|l . we obtain four equations for the amplitudes Ai^(Kk) (i = 0, 1, 2, 3) 
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^{K n A 0n + n a A aa - iQAi a } = 0, (58) 
y^{/s n A ln + - iQ(-A 0n + A 0a )} = 0, (59) 

^{Kn^2 n + Ka^2a} = 0, (60) 
y^KnAsn + ^a^3a} = (fs/lb), (61) 

where the summation is performed over the eigenvalues n = ±, a — a, 6; that is, the indices n and a correspond to 
the short-range and long-range eigenfunctions, respectively. Solutions for these equations yield Eas. ((23(1 - 1)25(1 . 



VII. APPENDIX B. JOSEPHSON EFFECT 



We consider the limit k^L >> 1 when one can neglect the overlap of the short-range eigenfunctions corresponding 
to the eigenvalues k±. In this case solutions for Fq j i(x) may be represented in the form 

F (x) = ^(A „f 2 ± Aonh) exp(-K„(L =p x)) + A Qa T 2 cosh(K Q a;) + Ao a h sinh(ft Q x), (62) 

where n, a, as before, are equal to ± and a, b. The matrix F^lx) has a similar form. The amplitude of the long-range 
component of F 3 (x) is small. This can be seen from Eq. ((54(1 . The matrix F± (x) has a different spatial dependence 

Fi(x) = f 3 2j(A ln f 2 ± Ainh) exp(-K„(L =F %)) + A la r 2 cosh(K a x) + A\ a ri sinh(K a x), (63) 

The first term in Eas. ((62163(1 describes the modes fast decaying from the S/F interfaces at x — ±L and the second 
term corresponds to the LRTC. The coefficients Ai n ,Ai a and Ao n , Ao a are connected with each other by Eas. l(54l57|) 
(see Eas. 1(23(1 - (124(1 ). The additional terms Ao t i n , a Ti appear because the matrix fs has changed (see Ea. l(29() '). In order 
to find these amplitudes, one has to substitute the expressions Eas. 1(62163(1 into the boundary conditions 

df/dx + z(Q/2)f 3 [a u f] + \ x=±L = (fshb)(T2 cos(^/2) ± r x sin^/2)) • a 3 , (64) 
Performing the calculations in this way, we arrive at four equations for the matrices Fi(L) (compare with Eas. ((58(1 - 

CD) 

^^{tt n A 0n + ^aA)a sinh# a + iQAi a sinh# Q } = 0, (65) 
'^{{n n Ai n + K a A la cosh6» Q ) + iQ(A 0n + A Qa cosh6> Q )} = 0, (66) 

n , a 

^^{— z tsai'df^K n AQ n + i cot dn^A^ sinh# a } = 0, (67) 

^2{n n A 3n + K a A 3a cosh6» Q } = {fs/jb) cos(<^/2), (68) 

We expressed i*2n,a in terms of Fon,a making use Eas. ((23() - ((24() . The corresponding equations for the coefficients 
•^■i,n,a may be obtained in a similar way. These equations coincide with Eas. 1(65(1 - 1(68(1 if one makes the replacement 
sinh0 Q cosh#Q, and cos(^/2) — > sin(<y9/2). In the main approximation in the parameter Q/k/,. solutions for these 
equations are given by Eas. l(33|) - (|3Tjl . 
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VIII. APPENDIX C. QUASI-BALLISTIC CASE 

In this Section we represent formulas for the condensate function f(x) in the case of a strong ferromagnet or a 
large mean free path / when the condition (|41|l is fulfilled. Substituting Ea. <|53|) into Eq.JSTJ and performing the 
transformation (|llfl . we obtain the linearized Eilenberger equation for the new function f n (x) in the ferromagnet (for 
brevity we drop the subindex n) 

(iTsldf/dx + i(Q/2)[il [&!, f] + + a J - i(a h /2) [a 3 ,f] + = (/}, (69) 

where a u = 1 + 2|u|r, ah = 2/irsgnw. In order to solve this equation, we represent the matrix f(x) as a sum of 
matrices symmetric s(x) and antisymmetric a(x) in the momentum space 

f(x)=s(x) + a(x), (70) 

Equations for these matrices can be obtained if we substitute Eq.{7UJ) into Ea. (|fi5|> and split it into the symmetric 
and antisymmetric in /i parts. We write down, for example, equations for the diagonal elements of matrices s(x) and 
a(x) in the spin space: s(x)ii(22) = s±, a(x)n(22) = a±. These equation have the form 

{K% - i M 2 l 2 d 2 /dx 2 )s ± + {Qln) 2 (K±/a u )s = i{^l) 2 Q[l + K ± /a u ]hds x /dx + K±(s±), (71) 



K±a± = - {fj,l)[i 3 d x s± + iQh], (72) 

The coefficients K± are defined in Ea. (|46|l . We solve equations for s± assuming that the coefficients K± are large, 
that is, ah is large. In addition, we assumed that the ratio Ql/ah — Qv/h is small. In this case the second term 
on the left-hand side and all terms on the right-hand side of Eq.J7IJ can be neglected. Solving this equation with 
the boundary conditions H43|) . in the main approximation we obtain the expression (|46|l . One can easily check that 
the term (s±) is much smaller than K±s± provided the quantity ah is large. The equation for the matrix §i can be 
readily obtained in a similar way. We get 

(al + (^) 2 [Q 2 (^i) 2 - d 2 /dx 2 })(T 3 h) ~ (lfi) 2 iQd x {s (l + { — f) + i—h] = atufoh) (73) 
ah ah a h 

The boundary condition for f 3 si(x) at x — requires that a\ — at x = 0. The expression for the antisymmetric 
matrix a\ has the form 

a^ai = -(ln)(d(T 3 si)/dx + iQs ) (74) 

As follows from Ea. l46f) . at the S/F interface in the main approximation so(0) = 0. Therefore the boundary condition 
for the matrix t 3 si(x) may be written as 

cUi = 0. (75) 

Eq. I|73|) can be solved in the following way. In the main approximation the coordinate dependence of the matrices 
SoM x ) 1S gi yen by Ea. fiBT (,s ,3 = (s + ± s_)/2) and these functions vary over distances \l/K±\ » v/h that are much 
shorter than a characteristic scale for the LRTC variation. Therefore approximately we can represent a solution for 
Eq.(|73Jl in the form 

§x(x) = sih{x) + 6sx(x) (76) 

where si/i(a;) is a short-range part of si(x) which is determined by: (d / 'dx)sih(x) = iQT 3 [§o(x) + i(a UJ /ah)s 3 (x)]lii 
particular the matrix 5si(x) contains the LRTC. As follows from Ea. (|74|l and Ea. (|75|l . the boundary condition for 
the function dsi(x) is 
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d x 8si{x) = -Q[{olu / ah)hh{x)] x =Q (77) 

where s^ifi) = T^t^fs- The equation for 5s\(k) in the Fourier representation (S§i(k) — J dx5si{x) exp(ifcx)) can be 
easily obtained from Eas. l|73l7t)l77|l . It has the form 

(al + Q 2 (^) 2 + (M) a )«i(*) = ^(Shik)) + 2(lfx) 2 Q^r 3 s 3 (0) (78) 

Oih Oth 

From this equation one can easily find 5si(k) 

ss ^ = « ^iSm <4^ 3(0)) + »' hm (7S) 

where N(k,fi) = o? u + Q 2 {a bJ /a h ) 2 + {kfil) 3 and N LR (k) = 1 - a u {l/N(k, n)}. The behavior of the LRTC h(x) 
is determined by poles of the functions N(k,n) and N^R^k). The first function has poles at k w These poles 
determine a variation of the matrix §\(x) over distances of the order of the mean free path from the S/F interface. 
The function N^n(k) has poles at much smaller wave vectors k which determine a long-range penetration of the triplet 
component. Indeed for k « r 1 we have: N LR (k) = (l 2 /3)[k 2 + K^] with K 2 Q = 2\u\t/D + (Q/a h ) 2 . Therefore at 
these wave vectors k the expression for 5s\{k) is reduced to Eq.ljlTS)!. 
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